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Un indicador es un atributo cualitativo o cuantitativo que 
permite describir y evaluar un aspecto específico de un 
sistema. Debe ser sencillo, construible con datos accesibles 
y de fácil interpretación. Diferentes organismos y programas 
internacionales han logrado una compilación de indicadores 
que aportan información sintética del estado y evolución 
de los sistemas hídricos. Sin embargo, es posible que sea 
necesario recurrir a una adaptación de los mismos, a fin 
de poder reflejar las particularidades de cada ambiente. 
En este trabajo se presentan los indicadores definidos para 
monitorear la calidad del agua subterránea, única fuente de 
provisión de la ciudad de Esperanza (Santa Fe, Argentina). 
La ruptura del equilibrio hidráulico del sistema acuífero 
multicapa explotado para abastecimiento ha producido 
variaciones químicas significativas en el agua de suministro. 
La aplicación de los métodos de agrupamiento (jerárquico 
y por K-medias) y componentes principales ha permitido 
verificar la presencia de este proceso, determinar cuál es 
el grado de afectación y ocurrencia temporal en cada uno 
de los pozos. Los resultados muestran que, si bien existen 
modificaciones en los valores de conductividad eléctrica 
del agua (CE) —normalmente utilizado como indicador de 
calidad—, este parámetro no es lo suficientemente indicativo 
del deterioro de la calidad del agua del sistema acuífero que 
se explota. Otros atributos dan indicios del ascenso de agua 
salada subyacente y de flujos descendentes desde el acuífero 
libre suprayacente, respectivamente, por lo que se establecen 
valores umbrales y niveles de alerta para la CE, el residuo 
seco, cloruros, sulfatos, nitratos y calcio. 


Palabras clave: indicadores de calidad del agua subterránea, 
protección de pozos, análisis de agrupamiento. 
Introducción y objetivos 


Las aguas subterráneas constituyen una im- 
portante fuente de abastecimiento tanto en áreas 


Abstract 


Paris, M., Tujchneider, O., Pérez, M., € D'Elía, M. (July-August, 
2014). Supply Well Protection: Groundwater Quality Indicators. 
Water Technology and Sciences (in Spanish), 5(4), 5-22. 


An indicator is a quantitative or qualitative attribute used to describe 
and evaluate a specific aspect of a system. It must be simple, buildable 
with available data and easily represented. Different organizations 
and programs have developed a compilation ofindicators that provide 
synthetic information of the state and evolution of water systems. 
Nevertheless, these indicators may need to be adapted in order to 
reflect the particularities of each environment. This paper presents 
the indicators defined to monitor groundwater quality for the aquifer 
underlying the city of Esperanza (Santa Fe, Argentina), the only 
source of supply for the region. The disruption of the hydraulic 
equilibrium of the multilayer aquifer system that is being exploited 
has produced significant chemical changes in the water supply. The 
results show that, although the electrical conductivity (EC) values 
(usually used as a quality indicator) of the water have changed, this 
parameter is not a sufficient indicative of the deterioration of water 
quality in the exploited aquifer system. Cluster (hierarchical and 
K-means) and principal methods have made it possible to verify this 
process, to determine the affectation level and its temporal occurrence 
for each exploitation well. Other attributes show evidence of the 
rise in the overlying salt water and of downward flows from the 
overlying unconfined aquifer, resulting in values for EC threshold 
and alert levels, total dissolved solid, chlorides, sulfates, nitrates and 
calcium. 


Keywords: Cluster analysis, groundwater quality indicators, well 
head protection. 


urbanas como rurales. Sin embargo, la ausencia 
de un apropiado diagnóstico y/o pautas de 
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manejo dentro de un plan de gestión de los 
recursos hídricos conducen en más de una 
oportunidad al deterioro de los mismos. 

La necesidad de disponer de agua segu- 
ra (en cantidad y calidad) pone en evidencia 
la importancia de llevar adelante políticas y 
acciones que garanticen la protección de las 
aguas subterráneas. Para una administración 
ambientalmente confiable de las aguas subte- 
rráneas, la mejor práctica es proteger el recur- 
so antes de su contaminación. Algunos países 
tienden a optar por la creación de programas 
de protección de las aguas subterráneas ante la 
contaminación, basados en el establecimiento 
de perímetros de protección de pozos (PPP) 
y /o en la cartografía de la vulnerabilidad de 
los acuíferos (Hirata €£ Reboucas, 1999). En 
ambos casos, la efectiva protección de los po- 
zos de abastecimiento requiere considerar, en 
primer lugar, la complejidad hidrogeológica de 
la zona de estudio, debida especialmente a las 
relaciones hidráulicas y a las características hi- 
drodinámicas e hidroquímicas de los cuerpos 
acuíferos que lo integran. La correcta concep- 
tualización del acuífero aportará bases sólidas 
para: 


e Establecer distintos niveles de restricciones 
y controles a los uso del territorio, sean 
estos pasados, presentes o futuros, con una 
incertidumbre aceptable. 

e  Instrumentar las herramientas legales, eco- 
nómicas, administrativas y los cambios ins- 
titucionales necesarios, para implementar 
prácticamente las restricciones y los con- 
troles al uso del territorio. 

e Definir de manera correcta los esquemas de 
explotación: ubicación y diseño de las obras 
de captación, caudales y cronogramas de 
extracción, etc. Las áreas de protección 
dependen del caudal de explotación, 
en especial en acuíferos no confinados. 
Bombeos no planificados de modo correcto 
pueden producir la ruptura del equilibrio 
hidrodinámico en el sistema acuífero, con 
la consecuente afluencia de aguas de otra 
calidad. 


e Diseñar un esquema de monitoreo del 
sistema en explotación a través de medi- 
ciones sistemáticas de variables indica- 
doras, que reflejen las eventuales varia- 
ciones de la calidad y cantidad del recurso. 
Esto brindará una alerta temprana que 
permitirá encarar medidas correctivas y/o 
paliativas ante su deterioro. 


En el caso de la ciudad de Esperanza, la 
creciente demanda de agua, el peligro de 
contaminación de las perforaciones ubicadas 
en el área urbana y el deterioro de la calidad del 
agua subterránea por ascenso de agua salina 
fueron importantes motivos para reubicar los 
pozos de abastecimiento de la ciudad en un 
nuevo campo de bombeo emplazado en la 
zona rural, diseñado a partir de un modelo de 
gestión desarrollado en 1999, como producto 
de las investigaciones antecedentes de este 
equipo de trabajo (Filí, Tujchneider, Paris, 
Pérez, € D'Elia, 1999). Si bien la Operación 
de estos pozos se materializó a partir del año 
2004, ningún esquema de protección ha sido 
aún formalmente definido y llevado adelante, 
en especial si se tiene en cuenta la intensa 
actividad agrícola-ganadera e industrial que 
se desarrolla en la región, las posibilidades de 
flujos descendentes y/o ascendentes a través 
del discontinuo estrato acuitardo que sobreyace 
al acuífero semiconfinado en explotación, la 
ruptura del equilibrio hidráulico del cuerpo 
acuífero inferior de alta salinidad producida 
por el bombeo no controlado, así como el 
posible efecto de la importación de agua 
superficial, que en un futuro complementaría 
el abastecimiento. Por ello, el objetivo 
general de esta investigación ha sido definir 
indicadores que colaboren con el monitoreo y 
control de la calidad del agua subterránea en 
el área de estudio, partiendo de la hipótesis de 
que la medición de la conductividad eléctrica 
del agua no es suficiente para evidenciar de 
forma temprana el deterioro de estos recursos 
hídricos. 

Diferentes 


organismos, asociaciones y 


programas internacionales han logrado una 
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compilación y clasificación de diferentes 
indicadores (Aureli € Taniguchi, 2008; Berger 
€ lams, 1996; Foster €: Loucks, 2006; Vrba 
€ Lipponen, 2007; Webb, Hirata, Kruse, éz 
Vrba, 2006). En Argentina, Blarasin, Cabrera, 
Matteoda, Damilano y Giuliano-Albo (2008) 
presentan una aplicación para evaluar cambios 
ambientales en acuíferos; Bocanegra, Massone, 
Irutig, Cionchi y Panizos (2009), una sobre 
los esquemas de explotación. Sin embargo, 
es posible que sea necesario recurrir a una 
adaptación de los mismos, a fin de poder 
reflejar las particularidades de cada ambiente. 
En relación con la calidad del agua subterránea, 
debe tenerse especialmente en cuenta que 
su calidad es la resultante de la interacción 
del agua de recarga, de las reacciones con los 
sedimentos y del tiempo de contacto o tiempo 
de residencia en el sistema (Doménico, 1972). 
El agua subterránea es un agente geológico 
de carácter general. Esto se debe de manera 
fundamental a la interacción agua-geología, 
que involucra una serie de procesos químicos, 
físicos y cinéticos dados por una evolución 
natural que tiende a alcanzar un estado 
de equilibrio o de mínima energía y, a la 
característica jerárquica y sistematizada de las 
líneas de corriente, que definen al flujo del agua 
subterránea como mecanismo de transporte 
y distribución de los fenómenos producidos 
(Tóth, 2001). Según Vrba y Lipponen (2007), un 
indicador de aguas subterráneas es un atributo 
que permite evaluar y describir un aspecto 
específico de este sistema hídrico. Deben 
permitir de una manera sencilla y rápida 
realizar el seguimiento del estado de las aguas 
subterráneas o de las modificaciones que se 
produzcan (Pernia-Llera, Lambán-Jiménez, é 
Molinero-García, 2005). Puede ser cuantitativo 
o cualitativo, según las necesidades que deba 
cubrir, pero fundamentalmente sencillo, pues 
no sólo debe ser fácil de construir a partir de 
datos medibles u observables, sino también 
interpretable; constituyen un soporte a la 
gestión sostenible de los recursos hídricos, 
sobre el 
estado actual y evolución del sistema hídrico 


proveen información sintética 


subterráneo, colaboran en el análisis de los 
procesos naturales e impactos de la actividad 
antrópica —en tiempo y espacio—, y facilitan 
la comunicación y participación pública en 
el proceso de planificación y definición de 
políticas. 


Principales características del área de 
estudio 


La ciudad de Esperanza (Departamento Las 
Colonias) se encuentra en el centro oeste de la 
Provincia de Santa Fe (Argentina), a 31? 27' de 
latitud sur y 60% 55'de longitud oeste, con una 
altura promedio sobre el nivel del mar de 36 m y 
una población de 35 000 habitantes. El principal 
acuífero en explotación tiene buen rendimiento 
y calidad. Desde esta área se complementa el 
abastecimiento de agua a una localidad vecina 
—Rafaela, Departamento Castellanos—, 
mediante un acueducto de 50 km, puesto en 
funcionamiento a fines de la década de los años 
setenta (figura 1). Más de 250 000 personas 
se abastecen de agua subterránea por red 
domiciliaria o extracciones particulares en estos 
dos departamentos provinciales. Esto deriva en 
un volumen de explotación estimado superior 
a los 51 millones de litros por día, de los cuales 
más de 26 millones de litros se consumen sólo 
en las localidades de Esperanza y Rafaela. El 
principal sustento económico del sistema social 
delaregiónlo brinda la agricultura, la ganadería 
que 
acciona de modo permanente sobre el sistema 


y una variada actividad industrial, 


acuífero, modificando sus características cuali- 
cuantitativas (Paris, Tujchneider, Pérez, D'Elia, 
éz Fili, 1998). En la actualidad, no se cuenta con 
un registro de regantes que permita estimar 
el volumen de agua subterránea utilizado 
para riego de los cultivos. Las explotaciones 
agropecuarias utilizan el agua subterránea para 
consumo humano, pecuario y riego de cultivos 
(D'Elia, Tujchneider, Paris, Pérez, € Pusineri, 
2012). No obstante, con datos del Censo 
Nacional Agropecuario 2002 (realizado por el 
Instituto Nacional de Estadísticas y Censos), se 
ha podido establecer que en el Departamento 


Tecnología y 


Ciencias del Agua, vol. V, núm. 4, julio-agosto de 2014 


Ciencias del Agua, vol. V, núm. 4, julio-agosto de 2014 


«o180j0u3aJ 


Paris et al., Protección de pozos de abastecimiento. Indicadores de la calidad del agua subterránea 


Argentina 


TF 
| ! 


e] Brasil] 


—— 


pl = 
$0 


ma Pacífi. 


Océa 


/ 


1 0 
HA 


Buenos Aires 


€ Pozos de explotación 
e Pozos de monitoreo 


1 kilómetros 


Figura 1. Ubicación del área de estudio. Localización de los pozos para abastecimiento de agua. 


Las Colonias se practica riego con agua 
subterránea mediante la explotación de 13 
perforaciones. Los 10 establecimientos rurales 
que realizan dicha práctica cultural están 
dedicados principalmente a la producción de 
maíz, soja y hortalizas. Aunque no se dispone 
de otra información más actualizada, se estima 
que en función de la producción agrícola y los 
periodos hídricos deficitarios registrados, a la 
fecha, la cantidad de regantes es mayor y, por 
ende, su porcentaje de utilización más elevado. 

El servicio de agua potable en la ciudad de 
Esperanza data del año 1930, cuando Obras 
Sanitarias de la Nación (OSN) construyó 
una red de distribución alimentada por tres 
pozos que producían alrededor de 180 m*/h 
(4 320 m*/día). Posteriormente, en 1940, se 
incorporaron dos pozos más; en 1969, otro; 
en 1974 se sumaron dos, y nueve en 1988. En 
1999 se contabilizaban 18 perforaciones en el 
área urbana, de las que solamente ocho estaban 
en funcionamiento, como consecuencia de 
problemas vinculados con el deterioro de la 


calidad del recurso (aumento de los contenidos 
de cloruros, salinidad y/o nitratos) (Paris, 
Tujchneider, D'Elia, Pérez, € Fili, 1999; Paris, 
2005). 

Hacia fines del año 1994, el servicio 
sanitario de agua y cloacas en las ciudades de 
Esperanza y Rafaela (como en otras localidades 
de la provincia) estuvo a cargo de una empresa 
privada concesionaria (Aguas Provinciales 
de Santa Fe Sociedad Anónima, APSF). Hasta 
aproximadamente el año 1997, la empresa 
concesionaria continuó con el esquema de 
explotación anteriormente descrito. Pero la 
creciente demanda, los peligros de contami- 
nación asociados con la actividad económica 
del área (localizada en especial en el sector 
urbano) y la falta de obras de saneamiento, 
condujeron al continuo deterioro de la fuente de 
abastecimiento (Tujchneider, Paris, Fili, D'Elia, 
éz Pérez, 1998; Paris et al., 1998). Desde 2006, 
el servicio sanitario (agua potable y cloacas) 
es suministrado por una empresa privada, 
donde el estado provincial es mayoritario en 
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su propiedad (Aguas Santafesinas Sociedad 
Anónima, ASSA). 

Dada la que 
abastecimiento de agua potable a estos dos 


importancia reviste el 
importantes centros urbanos a través de una 
fuente de origen subterráneo, su caracterización 
es fundamental en relación con la adecuada 
gestión del recurso en una región. Por ello, 
las investigaciones se han focalizado en el 
análisis del funcionamiento de este sistema y 
en el comportamiento de la calidad del agua 
subterránea en las áreas de mayor explotación 
histórica del acuífero: ejido urbano de la 
ciudad de Esperanza y campo de bombeo en 
el área rural. La figura 1 muestra la ubicación 
de las perforaciones que hoy día conforman 
el sistema de abastecimiento de ambas 
localidades (algunas de ellas actualmente 
fuera de servicio), pozos de estudio y pozos 
de monitoreo. Este nuevo campo de bombeo 
unificado está totalmente emplazado en el área 
rural del distrito Esperanza. 

Con base en conceptualizaciones antece- 
dentes (Filí et al., 1999; Tujchneider et al., 1998; 
Tujchneider, Pérez, Paris, £z D'Elia, 2005) y la 
interpretación de perfiles litológicos y estudios 
disponibles, Paris (2010) elaboró el esquema 
del modelo conceptual de funcionamiento del 
sistema hidrogeológico en el área de estudio 
que se presenta en la figura 2. El estrato 
superior, con un espesor variable entre 15 y 25 
m, corresponde al Grupo Pampa (sedimentos 
pampeanos) (Pleistoceno  lacustre-palustre- 
eólico), definido por Tujchneider (2000). Se 
trata de limos, arcillas y loess, pardos claros, 
oscuros, rojizos o grises, según el ambiente 
de depositación, con presencia de tosca en 
algunos niveles. Aloja un acuífero libre de bajo 
rendimiento y calidad variable. En la porción 
más profunda, donde los sedimentos son más 
arcillosos (de 2 a 3 m de espesor promedio), 
se comportan como un manto acuitardo 
discontinuo (Filí et al., 1999; Tujchneider et 
1998; Dalla-Costa, Díaz, 
2007). Inferiormente se localiza la formación 


al., € Aceñolaza, 


Ituzaingó (“arenas puelches”) (Plioceno 


fluvial), compuesta por arenas finas, medianas 


y gruesas amarillas, con recubrimiento 
limonítico. Espesor: 25-35 m; T = 600 m?/ 
día. Aloja un acuífero semiconfinado de buen 
rendimiento y calidad. Subyaciendo a ellas, la 
formación Paraná (Mioceno marino): porción 
cuspidal, arenas grises, arcillas arenosas y 
arcillas verdes, agua de elevada salinidad, 
a una profundidad de entre 45 y 53 m. Para 
Tujchneider et al. (2005), el comportamiento 
hidráulico de este sistema es de tipo multicapa, 
con la posibilidad de flujo descendente 
y/ O ascendente a través del acuitardo, en 
función de las relaciones de carga hidráulica 
imperantes. La extracción de agua del acuífero 
semiconfinado puede inducir no sólo el flujo 
desde el acuífero libre suprayacente, sino 
también incrementar la afluencia lateral de 
agua del propio acuífero e inducir el ascenso de 
agua salada. Como consecuencia del bombeo 
no controlado, la superficie piezométrica acusa 
en el año 1994 considerables descensos en el 
sector este (Tujchneider et al., 1998; Paris et al., 
1999). Dos años más tarde, en 1996, estas líneas 
mostraron una considerable variación en su 
diseño, traducida en un descenso promedio 
de casi 2 m y una divisoria definida en la línea 
equipotencial de 29 m. Hacia 1999 se logra una 
recuperación paulatina de los niveles, producto 
del cambio en el esquema de explotación. 

El sentido general de escurrimiento 
subterráneo es de oeste a este; la recarga local 
al primer nivel acuífero está en el orden de los 
63 mm/año y al segundo nivel, a través del 
acuitardo, es de aproximadamente 18 mm. El 
análisis de 9'*O y 9%H en el agua subterránea en 
relación con la recta meteórica muestra que la 
recarga proviene de la precipitación local, que 
ocurre en los meses de verano y otoño, y no 
refleja procesos de evaporación significativos. 
Asimismo, los valores de tritio detectados 
indican que el tiempo de residencia del agua en 
el acuífero semiconfinado es mayor a 50 años y 
que en el acuífero libre resulta de una mezcla 
entre recarga submoderna y reciente. En 
particular, los contenidos de tritio detectados 
en uno de los pozos de monitoreo del acuífero 
semiconfinado sugieren la mezcla de agua 
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Figura 2. Esquema del modelo conceptual de funcionamiento del sistema hidrogeológico en el área de estudio. 


de los dos niveles acuíferos. La presencia de 
nitrato estaría corroborando esta situación, 
favorecida además por la discontinuidad 
del manto acuitardo (D'Elia, Tujchneider, 
Paris, 6 Pérez, 2007; D'“Elia, Tujchneider, 
Paris, Pérez, € Gervasio, 2008). De acuerdo 
con contribuciones antecedentes, las aguas 
son de tipo bicarbonatado sódico (Filí et al., 
1999; Tujchneider et al., 1998; Tujchneider et 
al., 2005). Sin embargo, desde mediados de la 
década de los años noventa, en algunos pozos 
de la zona urbana se evidenciaron variaciones 
que derivaron en el cambio de tipo de agua 
bicarbonatada hacia el tipo clorurada sódica. En 
particular, en la ciudad de Esperanza, durante 
el periodo 94-96, se pudo comprobar el efecto 
que la mayor extracción produjo en el descenso 
de niveles y en la salinidad del agua (expresada 
en términos de conductividad eléctrica). La 
figura 3 muestra los cambios acontecidos en la 
calidad del agua en un pozo que fue sometido 
a una intensa extracción de caudales; las 
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flechas señalan la modificación sustancial en 
los contenidos aniónicos, que derivan en el 
cambio de tipo de agua, según la clasificación 
de Piper-Hill. La figura 4 pone en evidencia 
el aumento de salinidad registrado en este 
punto de observación a lo largo del tiempo. La 
disminución de salinidad registrada hacia fines 
de 2003 es el resultado de la reprogramación de 
su explotación, articulada al nuevo modelo de 
gestión (Paris, 2010). 


Metodología 


La base de datos utilizada se conforma con 
437 registros completos de ocho variables 
físico-químicas (conductividad eléctrica CE, 
alcalinidad total TAC, cloruro CL, sulfato 
SO4?, nitrato NO3,, calcio Ca”, magnesio 
Mg”, sodio+potasio Na'+K*). Estos registros 
corresponden a determinaciones analíticas 
practicadas en muestras de agua extraídas 
durante el periodo febrero de 1990 a diciembre 
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Figura 3. Modificaciones en el tipo de agua, según la clasificación de Piper-Hill en el pozo el7. 
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Figura 4. Variación de las concentraciones de cloruro + sulfato versus sodio + potasio para el pozo el2. 
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de 2007 en: a) pozos de explotación de la ciudad 
de Esperanza y acueducto a la ciudad de 
Rafaela (de 1990 a 2004); b) pozos de monitoreo 
ubicados en el distrito Esperanza (desde 2002), 
y Cc) pozos de explotación del nuevo campo de 
bombeo unificado (a partir de 2004). Esta base 
de datos integra información proveniente de 
diversas fuentes: Organismos provinciales, 
empresas prestatarias del servicio, universidad, 
etcétera. 

El error admisible en el balance iónico fue 
del 10%, límite definido considerando los 
valores de salinidad de las muestras de agua. 

Se aplicaron Métodos Estadísticos Multi- 
variados (MEM) para identificar o verificar pro- 
cesos que caracterizan la calidad del agua (en 
función de las variables que cuantitativamente 
los expresan) (Furi, Razack, Abiye, Kebede, éz 
Legesse, 2012; Izady et al., 2012; Daughney, 
Moreau-Forunier, $  Raaij, 
2012), conocer cuál es el grado de afectación 


Morgenstern, 


de estos procesos en cada uno de los pozos 
y establecer sus periodos de ocurrencia (o 
variación temporal) (Paris, 2010). Para ello, 
las concentraciones son expresadas en meq 1, 
estandarizados, según la tipificación de la 
Ley Normal y normalizadas por cálculo 
de logaritmo natural (Paris, 2005). Así, este 
conjunto volumétrico de datos tiene como 
dimensiones: 


e Los pozos que hacen referencia a las 
muestras de agua obtenidas en puntos 
georreferenciados del espacio (a modo de 
objetos o unidades de observación). 

e Las variables físico-químicas. 

e El tiempo que hace referencia a la fecha en 
la que fueron realizados los muestreos. 


La figura 5a muestra una esquematización 
de la base de datos volumétrica y su 
representación algebraica. 

Se advierte que los periodos de tiempo 
de registro no son regulares y que no cuentan 
con igual cantidad de datos. Esto se debe a 
que no todos los pozos han sido muestreados 
de forma simultánea o porque no estaban en 


operación aún o habían salido del servicio. Con 
ello, el cubo de datos no es “sólido”, sino que 
presenta “oquedades”. Estas últimas adquieren 
considerable importancia, según el modo de 
procesamiento seleccionado. 

El procesamiento se realizó aplicando 
una transformación previa, que permite re- 
definir algebraicamente la base volumétrica 
anteriormente descrita en un “conjunto plano 
de datos”, que mantiene intrínsecamente las 
dimensiones originales (y, por lógica, toda la 
información) (figura 5b). 

Esta transformación o “aplanamiento” 
permite identificar O verificar procesos que 
caracterizan la calidad del agua (en función 
de las variables que cuantitativamente los 
expresan), conocer cuál es el grado de afec- 
tación de estos procesos en cada uno de los 
pozos (determinando con ello las áreas de 
ocurrencia, dado que los pozos son objetos 
espacialmente georreferenciados) y establecer 
sus periodos de ocurrencia (o variación 
temporal). En este sentido, se hace notar que, 
en rigor, un análisis completo del cubo de datos 
debería incluir la sucesiva parametrización de 
cada una de las dimensiones (Paris, 2010), lo 
cual puede resultar en una implementación e 
interpretación engorrosa y limitada. Por otra 
parte, puede suceder que la dimensión final de 
los bloques resultantes no sea la adecuada, en 
términos de significación física y lo algebraica, 
tal como sucede cuando en un periodo de 
tiempo se han muestreado pocos pozos o 
cuando los muestreos no se han realizado en 
idénticos periodos. Además, la selección de 
qué variable, objeto o tiempo fijar no es un 
problema menor. Por ejemplo, sería posible 
pensar que si se registró un deterioro de la 
calidad del agua por aumento de la salinidad, 
la CE podría ser la variable más indicativa de 
estos cambios, pero esta selección podría estar 
enmascarando otros resultados o despreciando 
otras variables significativas. 

Se utilizó el análisis de agrupamiento (AA) 
por aglomeración jerárquica en Modo R y Q. 
En cuanto a la medida utilizada para expresar 
las similitudes existentes en el conjunto de 
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Pozos 


Variables 


Aj Kk 


i=l.N 


j=1..M k=1...P 


a) Esquematización de la base volumétrica de datos, su representación algebraica 


B (ii, jj) 


ii=1...(NxP) ji=1...M 


b) Esquematización y representación algebraica de la transformación a conjunto plano 


Figura 5. Base de datos, volumétrica de datos y transformación por aplanamiento. 


datos experimentales, se indica que en el Modo 
R, la matriz de similitud fue cuantificada por 
el coeficiente de correlación y en el Modo Q, 
por el coeficiente de distancia euclideana. En 
ambos casos, las asociaciones son definidas 
por agrupamientos 
por promedio de pares 
ponderados como criterio de ligamiento. Para 


de tipo aglomerativo 
jerárquico, no 
la validación de los resultados logrados se 
utilizaron el AA no jerárquico por K-medias y 
el análisis de componentes principales (ACP) 
en Modo R y Q. El empleo de metodologías 
de procesamiento con distinto fundamento 
estadístico y matemático pretende disminuir 
la incertidumbre que puede traducirse a los 
resultados debido a falencias en la definición 
de la estructura de estos modelos (Paris, 2005). 


En el ACP Modo R también se partió de la 
consideración de la matriz de correlación. 
En cambio, en Modo O, las similitudes se 
cuantificaron con el coeficiente coseno de theta 
o índice de similitud proporcional, definido por 
Imbrie y Purdy (1962) (en Jóreskog, Klovan, éz 
Reyment, 1976). La retención de componentes 
se efectuó de acuerdo con los siguientes 
criterios: varianza total representada (superior 
al umbral del 70%); magnitud y representación 
de los autovalores (valor mayor a 1 y/o 
varianza representada individualmente mayor 
a 10%); valor de la comunalidad o porcentaje 
de reconstrucción (superior al 60%), y la 
relación autovalor versus orden, reflejada en el 
punto de quiebre en el diagrama scree. En el AA 
no jerárquico, las asociaciones fueron definidas 
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conforme a los resultados aportados por el AA 
jerárquicos y el ACP. Finalmente, con algunos 
por 
y/o componentes, se elaboran gráficos de 


registros categorizados agrupaciones 
Na*+K* vs. CI+S504* para complementar la 
interpretación de los resultados obtenidos. 
Como resultado de esta etapa se definen tres 
grupos que reflejan las consecuencias sobre la 
calidad del agua subterránea de los esquemas 
de explotación aplicados. Para estos conjuntos 
se calculan los estadísticos y se evalúa el ajuste 
a la ley de distribución normal de cada una de 
las series químicas. En función del análisis de 
esta información sintética se seleccionan los 
indicadores, sus valores umbrales y los niveles 
de alerta que señalan los cambios significativos 
en la calidad. 


Resultados y discusión 


El AA Modo R pone en evidencia la fuerte 
relación entre los valores de CE y CI, a pesar 
de que las aguas son predominantemente de 
tipo bicarbonatado (figura 6). Por otro lado, 
también se manifiesta una estrecha relación 
entre la salinidad (expresada por CE) con 
SO? y Na*+K*, lo cual se asocia con el carácter 


sódico del agua. Esta asociación expresa el 
conocido aumento de los iones de CI y SO? 
cuando las aguas son más salinas debido al 
Kps (producto de la solubilidad) de las sales 
de estos aniones, en especial el cloruro de 
sodio, cuya solubilidad es independiente de 
la temperatura. Este grupo se completa con el 
núcleo Ca?'-Mg?”* (relacionado a la dureza del 
agua) que conforma un subgrupo con NO), 
Si bien no se dispone de valores de pH que 
informen sobre las condiciones de acidez- 
alcalinidad del medio, la nitrificación traería 
aparejada una disminución del pH y ello favo- 
recería la disolución de carbonatos de calcio y 
magnesio, con el consecuente incremento de la 
dureza (Martínez, Bocanegra, $2 Costa, 1997). 
La TAC se muestra como el único elemento 
que conforma el grupo principal restante. El 
AA en Modo Q diferenció tres asociaciones 
principales, al considerar el nivel de corte 3.5. 
Por razones de escala, en la figura 6 se ilustra la 
cantidad de elementos (1) queintegra cada nivel 
de asociación (G). Los grupos definidos por 
esta metodología estadística de agrupamiento 
jerárquico se comparan con los obtenidos con 
ACP en Modo R. Las asociaciones logradas 
con este procedimiento son muy coincidentes. 


Cr CE SsO4 Naf+k' NO, Ca”"Mg" TAC 


Si 
ÉS 


GH G2 Y 


G4 G5 G6 G7 G8 G9 Gl0 Gll Gl2 G13 Gl4 G 


| 
/ 


Gl6 G17 G181 G19 (20 


a=1A n=6 120 118 n=2 117 n-49 n-7 n-3 n-26 n-4 n=4 n-15 n-11 475 n-2 n-10 n-3 Y n-1 f-2 


E 
gos 5 
3 3 3 
bil ¡9 
u 
a 3 
8 0.6 3) 
u 3 
2 g 2 
E 04 E 
5 A 
202 1 
D 0 
Grupo G1 


(a) Dendrograma modo R 


Grupo G3 
Grupo G2 


(b) Dendrogama modo Q 


Figura 6. Dendrogramas. 
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Se definieron tres componentes principales 
de acuerdo con los criterios de retención de 
componentes mencionados anteriormente. 

La figura 7 muestra los autovalores y la 
varianza por ellos representada en forma 
individual y acumulada y, las cargas que 
caracterizan a cada una de las componentes 
(Comp) en relación con las variables químicas 
consideradas. Se han sombreado con gris los 
valores de mayor peso. Esta configuración 
también se verifica en el gráfico de scores 
resultante del ACP en Modo OQ, que se muestra 
en la misma figura para las componentes 1 y 
2. Como puede observarse, la conformación 
de componentes arroja la misma asociación de 
variables definidas por el AA Modo R. 

En cuanto a la asociación de registros en 
los pozos, la figura 8 representa los registros 
que integran cada una de las tres componentes 
en forma comparativa con las asociaciones 
determinadas por los agrupamientos jerár- 


quicos y no jerárquicos (K-medias). En el caso 
del ACB los pozos en donde se evidencia 
afectación de una componente han sido 
identificados en función de los scores calculados 
en Modo R. Estas áreas de afectación se 
establecen considerando las marcas superiores 
al valor del desvío estándar de cada serie de 
scores (Lawrence € Upchurch, 1976). Este tipo 
de representación gráfica de los resultados del 
análisis de componentes principales en Modo 
R no es común, pues usualmente los scores se 
representan en mapas donde se distinguen las 
áreas de mayor afectación de las componentes 
como aquellas que exhiben valores superiores 
a 1 (en algunos casos se consideran todos los 
valores positivos). 

Con el razonamiento realizado respecto al 
“aplanamiento” del cubo de datos se debería 
generar un conjunto de mapas (uno para cada 
periodo de observación), cuya interpretación 
sería muy dificultosa. Debe tenerse en cuenta 


Autovalores y varianza representada - ACP Modo R ] 
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Componente 2 


Gráfico de scores - ACP Modo Q 
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Figura 7. Resultados ACP. 
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Registro en los pozos 


Figura 8. Agrupamientos logrados por métodos jerárquicos, no jerárquicos (K-medias) y análisis de componentes 


principales (Modo R). 


nuevamente que no todos los pozos cuentan 
con muestras en todos los periodos, por lo 
que la densidad de información en cada 
mapa sería muy variable y su evaluación 
muy limitada. Además, no sería factible hacer 
una comparación con los resultados de los 
agrupamientos anteriormente presentados. 

Se desprende de la figura 8 que algunos 
registros en los pozos no tienen representati- 
vidad en ninguna de las tres componentes. 
Esto se debe a las propias características del 
análisis de componentes principales, donde las 
componentes se seleccionan garantizando un 
nivel umbral de varianza representada. Esta 
es una diferencia sustancial con las metodolo- 
gías de agrupamientos que condicionan la 
asociación de elementos al nivel de corte elegido 
(método jerárquico) o a la cantidad de grupos 
predeterminada (K-medias). Sin embargo, 
los resultados son consistentes, porque los 


registros de los pozos más representativos de 
cada uno de los grupos son coincidentes. 
Luego de realizar esta validación cruzada de 
los resultados de diferentes MEMs, en el cuadro 
1 se presentan los valores medios y medidas 
de la tendencia central para las variables 
consideradas, según los grupos principales 
definidos. El grupo G1 engloba las muestras 
de agua correspondientes a los pozos donde se 
han seguido cronogramas de explotación que 
derivaron en valores elevados de salinidad. 
En este conjunto, la variación temporal de Ct 
+SO,? es simultánea a la de Na'+K*. El grupo 
G2 comprende los pozos de explotación y de 
monitoreo en periodos donde se manifiesta 
un aumento paulatino de la salinidad, aunque 
sin ser notorio en la expresión de la CE. En 
este caso, los contenidos salinos se modifican 
a lo largo del tiempo en términos de SO? y 
Na'+K*, pero manteniendo un umbral de 
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CI. Se observa también un aumento de las 
concentraciones de Ca”. Una causa posible 
del aumento de los iones mencionados puede 
atribuirse al ascenso de agua clorurada sódica 
alojada en los sedimentos de origen marino 
de la formación Paraná, al aporte de sulfato 
de calcio proveniente de los cristales de yeso 
y arcillas arenosas calcáreas de la misma 
formación (Iriondo, 2007; Tujchneider, 2000) 
en la búsqueda del equilibrio ante el cambio en 
las concentraciones (y posiblemente presiones) 
que experimenta el sistema (Principio de Le 
Chátelier) (Appelo é: Postma, 1993). En el 
grupo G2, las concentraciones de NO, y Mg” 
son superiores a las del grupo G3, pero menores 
a los estadísticos del grupo Gl. Esto podría 
estar indicando el aporte de agua alojada en 
el acuífero libre sobreyacente, especialmente 
en las áreas donde el manto acuitardo podría 
presentar discontinuidades. A su vez, como se 
indicó antes, la nitrificación traería aparejada 
una disminución del pH y ello favorecería la 
disolución de carbonatos de calcio y magnesio, 
con el consecuente incremento de la dureza 
(Martínez et al., 1997). Por su parte, el grupo 
G3 refleja el conjunto de muestras de agua no 
afectadas aún por el proceso de salinización 
(aumento de sales de cloruros y sulfatos); 
los valores de CE no muestran aumentos 
importantes en relación con los registrados 
al inicio de la operación de los pozos y los 
cambios temporales en los contenidos iónicos 
no son importantes. Es decir, estos grupos se 
diferencian no sólo por el aumento de la CE (G1 
vs. G2 y G3) sino también por el incremento de 
los contenidos de CI, SO,?, NO,, Ca” y Mg”. 
Por otra parte, se señala que los tres grupos 
identificados por el procesamiento estadístico 
multivariado de las series química reflejan las 
consecuencias de los esquemas de explotación 
a los que fuera sometido el acuífero. Así, por 
ejemplo, en función del caudal bombeado, una 
misma perforación puede contar con registros 
en los distintos grupos. Tal es el caso que ilustra 
la figura 9 para el pozo el7. Según este gráfico, 
los estadios temporales comprenderían los 
periodos 1991-1995, 1995-1999 y 1999-2006. Si 


bien no se dispone de registros de caudal para 
este lapso, se señala que tal comportamiento 
químico es notablemente concordante con 
las variaciones sufridas por la superficie 
piezométrica anteriormente descritas. Así, la 
aplicación de los MEM en la identificación y 
verificación del proceso de ascenso de agua 
salada ha determinado el grado de afectación 
en cada uno de los pozos y su ocurrencia 
temporal. 

De acuerdo con ello, los grupos Gl, 
G2 y G3 engloban las muestras de agua 
correspondientes a los pozos, donde: a) se 
han seguido cronogramas de explotación que 
derivaron en el incremento de la salinidad 
(salinizados); b) se manifiesta un aumento 
paulatino de la salinidad (estadio intermedio), 
y c) pozos no afectados aún por el proceso de 
salinización (línea base), respectivamente. 

A partir de aquí se plantea la definición de 
indicadores que permitan poner en evidencia 
en forma temprana los cambios en la calidad 
del agua anteriormente descritos. Así, la 
definición del indicador deberá considerar 
los valores 
característicos de la composición natural del 


el fondo natural, es decir, 
agua a partir del cual se pondrá en evidencia 
el estado del sistema, las presiones a que esté 
sometido y su respuesta. Manzano, Custodio 
y Nieto (2003) definen al fondo natural de la 
calidad del agua subterránea como “el rango 
de concentraciones de un elemento, especie o 
sustancia dados, presentes en fase disuelta y 
derivado de fuentes geológicas, biológicas o 
atmosféricas naturales”. Los mismos autores 
señalan además que el proyecto BaSeLiNe 
(Natural BaSeLiNe Quality in European 
Aquifers) decidió adoptar operativamente a 
la mediana como parámetro más ilustrativo 
del valor característico del fondo natural del 
agua subterránea y los percentiles de 2.3% y 
97.7% para ilustrar su rango de variación. Esto 
asegura que el 95.4% de la población estudiada 
esté dentro del rango. Sin embargo, los 
mismos autores señalan que si la distribución 
estadística es normal o lognormal, la media o 
la mediana podrán considerarse como valores 
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Figura 9. Evolución de CE y SO,” en el pozo el7 y su relación con las asociaciones determinadas. 


característicos, pero con frecuencia esto no 
ocurre debido a que la calidad del agua es el 
resultado de más de un proceso. Investigaciones 
realizadas por Blarasin, Damilano, Cabrera y 
Matteoda (2006) indican que la distribución 
fl presenta un mejor ajuste a la distribución 
empírica y determinaron un estimador m* del 
fondo natural, que resulta de una combinación 
lineal de la media y la mediana. Señalan, 
no obstante, que es necesario estudiar con 
mayor profundidad las distribuciones que 
siguen las distintas componentes si el objetivo 
es determinar el fondo natural a través de la 
estimación paramétrica de las concentraciones 
típicas (parámetros de localización) de cada 
componente y su distribución en torno a estos 
valores. El análisis realizado en este acuífero 
sobre la base de los resultados presentados 
permite definir los indicadores de calidad del 
agua de abastecimiento que se muestran en el 
cuadro 1b. Para el caso de estudio, los valores 


de TAC son considerablemente uniformes (con 
un leve aumento a medida que disminuye 
la salinidad) por lo cual este parámetro 
no es útil para realizar una discriminación 
en el conjunto. Si bien se visualiza a la CE 
como una variable indicativa del ascenso de 
agua de mayor salinidad proveniente de la 
formación Paraná —que además es de fácil 
determinación, in situ y en tiempo real—, 
la interpretación del fenómeno de ruptura 
del equilibrio hidrodinámico se muestra con 
mayor elocuencia cuando se evalúa la relación 
CI+SO,?” vs. Na*+K*. Esto se manifiesta clara- 
mente en los pozos donde se comienzan a 
evidenciar cambios en la salinidad del agua 
(grupo G2), se modifican en términos de SO? 
y Na'+K*, pero mantiene los valores de CI. 
Por otra parte, la asociación entre Ca”, Mg? 
y NO, (este último también vinculado con 
TAC), manifiesta tanto en el AA como en el 
ACP (Modo R), y puede estar dando indicios 
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Cuadro 1. Construcción de indicadores de calidad del agua de abastecimiento. 


a) Valores medios y medidas de dispersión de las series químicas según los grupos 


CE ho [ell SO, NO, Ca? Mg” Na'+K* 
S/cm CO,Ca mg/l mg/l mg/l mg/l mg/l mg/l 
Grupo 1 Media 3172.22 515.67 427.56 549.83 19.67 41.61 24.50 682.76 
Desvío 615.29 30.07 WA 149.02 5.41 14.89 ES Asi 
Cuartil 1 2 535.00 497.50 347.50 406.25 17.25 29.00 17.75 578.03 
Mediana 3 390.00 515.00 490.00 605.00 21.00 45.00 27.00 747.13 
Cuartil 3 3 632.50 528.75 527.50 660.00 22.00 54.00 31.00 758.62 
Grupo 2 Media 1 388.63 542.79 84.96 94.98 33.26 24.76 So 297.93 
Desvío 280.77 48.63 50.26 52.75 32.27 12.37 6.28 70.87 
Cuartil 1 1225.00 515.00 38.50 50.00 8.75 17.50 10.00 268.75 
Mediana 1 399.00 541.00 90.50 100.00 23.50 24.00 13.00 299.50 
Cuartil 3 1540.50 570.00 108.88 135.00 42.00 30.25 16.00 341.96 
Grupo 3 Media 1363.25 549.81 69.93 103.99 10.47 13.08 7.82 311.16 
Desvío 161.23 897 24.82 34.81 8.69 4.92 2.69 38.88 
Cuartil 1 1252.00 520.00 56.00 81.00 6.00 10.00 6.00 288.00 
Mediana 1 368.00 545.00 68.00 100.00 8.00 12.00 7.00 317.00 
Cuartil 3 1 465.00 570.00 80.00 128.05 12.00 16.00 9.00 340.50 
b) Indicadores de calidad del agua de abastecimiento 
CE RSO Cr 10% NO, Ca? 
S/cm mg/l mg/l mg/l mg/l mg/l 
Umbral 1 368 2 062 68 100 8 12 
Alerta 1541 PIS85 109 135 42 30 


(*) Los valores de RS fueron estimados como en función de la función de regresión lineal entre los datos disponibles (RS = 1.5853*CE-106.97). 


c) Límites establecidos en la normativa de referencia para el abastecimiento de agua potable 


E) R 1 2 N - 2+ 

Ley 11.220 e S e SO, O, Ca 
S/cm mg/l mg/l mg/l mg/l mg/l 
Límite recomendado 1 000 250 100 25 100 
Límite obligatorio 1500 400 200 45 250 


(++) la ley provincial no fija límites para CE. 


de aportes del nivel acuífero superior (libre) 
y de cargas contaminantes provenientes del 
lixiviado de residuos ganaderos y fertilizantes, 
propias del área rural donde se sitúa el campo 
de bombeo. Esta situación se presenta con más 
elocuencia por los mayores valores de NO, 
del grupo G3, formado principalmente por 
muestras extraídas de pozos de explotación y 


de monitoreo localizados en el área rural. El 
valor umbral se define a partir de la mediana 
del grupo G3 y el nivel de alerta de cambios 
significativos con el tercer cuartil del grupo 
G2. El grupo G3 describe las condiciones de la 
línea base de la calidad del agua subterránea, 
considerando la información disponible. En 
cuanto a la CE, no puede dejar de reconocerse 
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su practicidad de medición en campo y en 
laboratorio, y que verdaderamente puede 
ofrecer indicios de salinización, pero como 
puede verse en el cuadro la, los tres cuartiles 
considerados no muestran diferencias en 
los grupos G3 y G2. Por lo tanto, no es un 
parámetro que ofrece una alerta temprana 
de manifestación del proceso ni tampoco del 
aporte desde el acuífero libre. Se incluyen los 
valores del residuo seco (RS), atendiendo el 
marco regulatorio para el abastecimiento de 
agua potable en la Provincia de Santa Fe (Ley 
Provincial 11220). El cuadro 1c muestra los 
límites de esta normativa. 


Consideraciones finales 


A fin de evitar la ruptura del equilibrio hidráu- 
lico del sistema acuífero multicapa que susten- 
ta el abastecimiento de dos importantes locali- 
dades de la provincia de Santa Fe (Argentina), 
debería realizarse el seguimiento de los indica- 
dores presentados en el cuadro 1b. Se trata de 
determinaciones químicas rutinarias, de fácil 
determinación, que permiten evidenciar las re- 
laciones entre el acuífero alojado en las “arenas 
puelches” (semiconfinado, donde se extrae el 
agua para consumo humano), el acuífero libre 
alojado en los sedimentos pampeanos y el agua 
salada alojada en la formación Paraná (de ori- 
gen marino). 

Se ha corroborado en esta investigación 
que, para el caso de estudio, la medición de la 
conductividad eléctrica del agua no es suficiente 
para detectar de modo temprano el deterioro de 
los recursos hídricos. Los indicadores definidos 
adquieren solidez y representatividad a partir 
del uso de diferentes métodos estadísticos 
multivariados (MEM) para la validación 
cruzada de los resultados y, fundamentalmente, 
del sustento que proporcionan otras evidencias 
concretas de la geometría y funcionamiento del 
sistema hidrogeológico. A su vez, los resultados 
logrados sirven de base para el planteo de líneas 
de investigación, que permitirán complementar 
el conocimiento de las relaciones hidráulicas 


entre los cuerpos acuíferos. Dado que el flujo 
vertical descendente, puesto de manifiesto por 
las variaciones de Ca”, Mg” y NO,, puede 
constituirse en el agente de transporte de 
contaminantes derivados del uso del territorio, 
en próximas etapas, mediante tareas de campo 
y gabinete, se evaluará el volumen de agua 
por unidad de área que está almacenado en 
esta capa y que puede ser liberado, teniendo 
en cuenta las características litológicas de estos 
acuíferos, el espesor de la capa acuitarda, su 
presencia, espesor y ausencia, y los parámetros 
(conductividad 


hidráulicos  formacionales 


hidráulica, transmisividad, coeficiente de 
almacenamiento, factor de goteo, etc.). En 
relación con el transporte de contaminantes 
hacia el ambiente subterráneo, se considerarán 
además las características del escurrimiento 
superficial en el área rural en las cercanías 
de pozos (por conducción en lámina y en 
cauce), y se evaluarán las prácticas de riego 
y particularidades del drenaje natural de 
los suelos. Por otra parte, considerando las 
variaciones detectadas en la calidad del agua 
debido al ascenso de agua desde el cuerpo 
acuífero inferior de alta salinidad producida 
por el bombeo no controlado y la consecuente 
ruptura del equilibrio hidráulico, se evaluará 
también por ensayos de bombeos, toma de 
muestras de agua, estudios geofísicos y otras 
tareas de campo, y por el análisis en gabinete 
de la información de caudales de extracción. 
Además se considera complementar las 
determinaciones químicas e isotópicas que 
permitan discriminar el origen de los aportes 


de nitrato en ambos cuerpos acuíferos. 
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